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We study the isospin effects of the mean field and two-body collision on 

the nucleon emissions at the intermediate energy heavy ion collisions by 

using an isospin dependent transport theory. The calculated results show 

that the nucleon emission number Nn depends sensitively the isospin effect 

of nucleon nucleon cross section and weakly on the isospin dependent mean 

field for neutron-poor system in higher beam energy region . In particular 

,the correlation between the medium correction of two-body collision and 

the momentum dependent interaction enhances the dependence of nucleon 

emission number Nn on the isospin effect of nucleon nucleon cross section. 

On the contrary, the ratio of the neutron proton ratio of the gas phase 

to the neutron proton ratio of the liquid phase ,i.e., the degree of isospin 

fractionation < {N/Z)gas >b / < {N/Z)iiq >b depends sensitively on the 

isospin dependent mean field and weakly on the isospin effect of two-body 

collision for neutron-rich system in the lower beam energy region. In this 

case, Nn and < {N/Z)gas >b / < {N/Z)iiq >b are the probes for extracting 

the information about the isospin dependent nucleon nucleon cross section in 

the medium and the isospin dependent mean field, respectively. 
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One can study the isospin effects of heavy ion colUsion to get the information 
about isospin dependent in-medium nucleon-nucleon cross section and symmetry 
potential. In order to obtain this information several interesting isospin effects in 
heavy ion collisions have been explored both experimentally and theoretically over 
the last few years[l — 16, 26 — 28]. For example, recently, Bao-An Li et al investi- 
gated the isospin effect of the mean field and showed that the neutron-proton ratio 
of preequilibrium nucleon emission,the neutron-proton differential flow and proton 
ellipse flow are the probes for extracting the isospin dependent mean fleld at lower 
beam energy region [1, 2, 14, 15, 16]. R.Pak et al have found that the isospin depen- 
dence of collective flow and balance energy is mainly originated from the isospin- 
dependent in-medium nucleon-nucleon cross section [17, 18]. However two essential 
ingredients in heavy ion collision dynamics, the symmetry potential and the isospin 
dependent in-medium nucleon-nucleon cross section have not been well determined 
so far. Recently Bao-An Li ,M. Colonna , M.Di Toro and V.Bar an et al studied the 
isospin fractionation in the intermediate energy heavy ion collisions in recent years 
[15, 16, 19, 30]. The isospin fractionation is an unequal partitioning for the neutron to 
proton ratio N/Z of unstable asymmetry nuclear matter between low and high den- 
sity regions. But the neutron-proton ratios in gas phase and liquid phase calculated 
by Bao An Li were from low densities less (free) and higher (bound)than pc—l/8po 
respectively because IBUU model can not calculate the formation of fragments[14]. 
Even though Bao An Li 's work proposed an useful physical point of review but it can 
not compare with the experimental data directly. An indications of this phenomenon 
have been found recently in intermediate energy heavy-ion experiments by H. Xu 
et al [9]. However they did not compare the calculated results with experimental 
data directly and there were some assumptions in their calculations for comparing 
with the experimental data. In this case, how to quantify the gas phase and the 
liquid phase in isospin fractionation process for comparing with the experimental 
data directly is important. Based on the isospin dependent quantum molecular 
dynamics model (IQMD)[2, 20, 21, 22, 26] we investigated the isospin effects of the 
nucleon emission number A^^ for neutron-poor systems in higher beam energy region 
and quantify the degree of isospin fractionation, (A^/Z)^^^ / {N/Z)iiq by using the 



neutron proton ratios of the gas phase and hquid phase for the neutron-rich systems 
in lower beam energy region. Where {N/Z)gas and {N/Z)iiq are the neutron-proton 
ratio of nucleon emission and that of the fragment emission respectively. The cal- 
culated results show that (l)the nucleon emission number depends sensitively 
on the isospin effect of nucleon nucleon cross section and weakly on the symme- 
try potential in the higher beam energies ,on the contrary, the degree of isospin 
fractionation {N/Z)gas / {N/Z)iiq is sensitive to the the symmetry potential and in- 
sensitive to the isospin effect of in-medium nucleon nucleon cross section. It is worth 
mentioning that our calculated results can be compared with the experimental data 
directly to get the information on the symmetry potential or the nucleon nucleon 
cross section in the medium. 

A quantum molecular dynamics (QMD) model [20,21,22] contains two dynamical 
ingredients,the density dependent mean field and the in-medium nucleon nucleon 
cross section. In order to describe the isospin dependence appropriately, the QMD 
model must be modified properly. The density dependent mean field should contain 
correct isospin terms, including the symmetry potential and the Coulomb poten- 
tial. The in-medium nucleon nucleon cross section should be different for neutron 
neutron(or proton proton ) and neutron proton collisions, in which Pauli blocking 
should be included by distinguishing neutrons and protons. In addition, the initial 
condition of the ground state of two colliding nuclei must also contain isospin effects. 
Considering the above ingredients, we have made important modifications in QMD 
to obtain an isospin dependent quantum molecular dynamics (IQMD) [2,26]. The 
initial density distributions of the colliding nuclei in IQMD are obtained from the 
calculations of the Skyrme-Hatree-Fock with parameter set SKM* [23], and the ini- 
tial code of IQMD is used to determine the ground state properties of the colliding 
nuclei, such as the binding energies and RMS radii, which agree with the experi- 
mental data for obtaining the parameters of interaction potential as an input data 
of the dynamical calculation by using IQMD. The interaction potential is given by 

, where f/^""' ^U'^^y , jjYuk^jjMDi jjPauh Coulomb potential, the density 

dependent Skyrme potential, the Yukawa potential, the momentum dependent in- 



teraction and the Pauli potential [31,32] respectively ( see in more detail in Refs 
[2,20,21,22]). 

In the present paper, three different forms of h.eive been used [1, 15], 

= cF^{u)St, (2) 
C/f = cF2{u)5t, + ^cF2{u)S' (3) 
Ur-cF,{u)[6r^-]s'] (4) 



C/o'^"^ = 0.0 (5) 



with 



{1 for neutron 
—1 for proton 

Here c is the strength of symmetry potential, taking the value of 32MeV, Fi(u)=u 
, F2(u)=u^ and F-i(u) = u^^'^,u= —,5 is the relative neutron excess 5 = ^"V^'' = 
^"p^^ -PiPoiPn Pp ^'^'^ the total , normal , neutron and proton densities, respec- 
tively. Where Uq'"^—Q.Q indicates no symmetry potential. It is worth mentioning 
that the recent studies of collective flow in heavy ion collisions at intermediate en- 
ergy have indicated a reduction of in-medium nucleon nucleon cross sections. An 
empirical density dependent nucleon nucleon cross section in medium [24] has been 
suggested as follows 

(7nn = (1 + a— )(7aw > (6) 
Po 

where (t^j^'^ is the experimental nucleon nucleon cross section [33]. The above ex- 
pression with a ~ —0.2 has been found to reproduce well the flow data. The free 
neutron proton cross section is about a factor of 3 times larger than the free proton 
proton or the free neutron neutron one below 400 MeV, which contributes the main 
isospin effect from nucleon-nucleon collisions. We used equation (6) to take into ac- 
count the medium effects, in which the neutron proton cross section is always larger 
than the neutron neutron or proton proton cross section in the medium at the beam 
energies in this paper. 



The isospin effect of the in-medium nucleon nucleon cross section on the 
observables is defined by the difference between the observables from an isospin 



dependent nucleon nucleon cross section ci"" and that from an isospin independent 
nucleon nucleon cross section a"'°^'^° in the medium. Here cr**° is defined as cr„p > 
Cnn—Cpp and cr""**" means (j„p = (7„„ = app, where (7„p, (T„„ and app are the neutron 
proton, neutron neutron and proton proton cross sections in medium respectively. 

In order to study the influence played by the medium correction of the isospin 
dependent nucleon nucleon cross section on the isospin effects of the nucleon emission 
, we investigated the number of nucleon emission Nn as a function of the beam 
energy for the mass symmetry neutron-poor system '^^Kr +'^^ Kj. (top panels) and 
mass asymmetry neutron-poor system ^^"^Sn Ca (bottom panels) with the same 
system mass at impact parameter b= 4.0 fm for the different symmetry potentials 
jj^ym ^ jjsym f/^*'™ as Well as the isospin dependent in-medium nucleon nucleon 
cross section a^'^" and the isospin independent one a"'°'^^" (the explanations about 
line symbols are in Figures) . It is clear to see that all of lines with flUed symbols 
are larger than those with open symbols,i.e.,all of Nn with a^^° are larger than those 
with cr""'*" because the collision number with u'*" is larger than that with cr""**". We 
also found that the gaps between lines with the filled symbols and the open symbols 
are larger but the variations among lines in each group are smaller. These mean 
that Nn depends sensitively on the isospin effect of in-medium nucleon nucleon cross 
section and weakly on the symmetry potential. In particular, the gaps between two 
group lines with a = —0.2 in left panels are larger than corresponding those with 
a = 0.0 in right panels, i.e., the medium correction of two-body collision enhances 
the dependence of Nn on isospin effect of two-body collision. 

To consider the contributions from all of impact parameters we also calculated 
an impact parameter average values of nucleon emission number < Nn >b for the 
same incident channel conditions and the same symbols as Fig.l. We can get the 
same conclusions as Fig.l. 

Why does the medium correction of nucleon nucleon cross section enhance the de- 
pendence of Nn on the isospin effect of two-body coUision? As well know that the 
isospin dependent in-medium nucleon nucleon cross section is a sensitive function of 
the nuclear density distribution and beam energy as shown in Eq.(6). Fig. 2 shows 
the time evolution of the ratio of nuclear density to normal one, — for the reac- 



tion ''^Kr + '^^Kr and symmetry potential f/-^^™ at E= 150 MeV/nucleon and b= 
4.0 fm. From the values of peak for ^ in the insert in Fig. 2 it is clear to see that 
p((7**°, a — —0.2) is larger than p{a'^'^^°, a — —0.2) and p{a^^°, a — 0.0) is larger than 
p(^^nmso^ Q, _ Q Q-^ because the larger coUision number from cr**" increases the nuclear 
stopping and dissipation , which enhances the nuclear density, compared to the case 
with or"°*'^°. From Fig.2 we can also see that ^'s decrease quickly with increasing 
the time after peak of ^. Because the larger compression produces quick expanding 
process of the colliding system ,while the small compression induces slow expanding 
process , at the same time, the ^'s decrease quickly with expanding process of sys- 
tem. Hovewor the decreasing velocity of ^ is larger for the quick expansion system 
than that for the slow expansion system, up to about 70 fm/c, on the contrary, 
p(^^nmso^ Q! = —0.2) is larger than p{a^'^", a — —0.2) and p{a^°^^°, a = 0.0) is larger 
than p{a'^^°, a — 0.0). In particular, the gap between two hues for a — —0.2 is larger 
than that for a — 0.0 after about 70 fm/c. This property is very similar to the N^, 
which means that the medium correction of an isospin dependent nucleon nucleon 
cross section enhances also the dependence of p on the isospin effect of two-body 
collision, which induces the same effects to iV„ through the nucleon nucleon cross 
section as a function of the nuclear density as shown in Eg. (6). 
We also found an important role of the MDI for enhancing the isospin effect of 
two-body collision correlating with the medium correction of two-body collision on 
the . Fig.3 shows the time evolution of Nn for a symmetry potentials f/*^'™ at 
beam energy of 100 MeV/nucleon and impact parameter of 4 fm. They are four 
cases (see Fig.3) for MDI in the left window and NOMDI in the right window. It 
is clear to see that the gap between the solid hue and dot hue with MDI in the 
left window is larger than corresponding gap with NOMDI in the right window in 
the medium(here only a = —0.2), i.e., MDI increases the isospin effect of two-body 
collision on the in the medium because above gap is produced from the isospin 
effect of nucleon-nucleon cross section in the medium. Here physically there are 
two mechanisms at work here. (1) The average momentum of a particle in medium 
is higher in a heavy ion collision than in cold nuclear matter at the same density. 
(2) MDI induces the transporting momentum more effectively from one part of the 



system to another, in which particles also move with a higher velocity for a given 
momentum than in free space. 

The degree of isospin fractionation [N / Z)gas/ {N / Z)iiq is measured by the ratio of 
the neutron-proton ratio of the nucleon emission (gas phase) to that of the interme- 
diate mass fragment emissions (liquid phase). Where the liquid phase includes all 

of fragments with the charge number from 2 to {Zp + Zt)/2 ,where Zp and Zt are 
the charges for projectile and target respectively 

Fig.4 is plotted the time evolution of the impact parameter average value of < 
{N/Z)gas >b / < {N/Z)uq >6(top window), < {N/Z)gas >b (in left side of bottom 
window) and < {N/Z)iig >b (in the right side of bottom window) for the different 
symmetry potentials U^^"^ and f/*^™ with the same a^'"^" or the same (t"°*'^° for the 
neutron- rich system ^'^'^Sn ^j^^ beam energy of 50 MeV/nucleon . It is 

clear to see that the variations for the < {N/Z)gas >b / < {N/Z)uq >6's with the (t'*" 
and a'^'^^° but the same symmetry potential U^^^ or JJi'^ are smaller, but the gap 
between the < {N/Z)gas >b / < {N/Z)iiq >6's with different symmetry potentials 
?7|^™ and f/-^^™ but the same (t"° or o-"°*'^° are larger. Namely, the amplitude of 

< {N/Z)gas >b / < {N/Z)iiq >6 depends sensitively on the symmetry potential and 
weakly on the isospin effect of in-medium nucleon nucleon cross section. 

It is also very clear to see that < {N/Z)gas >b has the same property as < {N/Z)gas >b 
I < {N/Z)iiq >b but < {N/Z)iiq >b is insensitive to both symmetry potential 
and nucleon-nucleon cross section in the medium. Namely above property for 

< {N/Z)n >b / < {N/Z)iiq >h is produced mainly from the gas phase . 

What is the main source for producing the isospin fractionation process in inter- 
mediate energy heavy ion collisions? As well know that in regions just off normal 
density the fields felt by neutrons and protons are very different for different symme- 
try potentials , especially, below normal density since symmetry potential is repulsive 
for neutron and attractive for proton. We thus expect important transport effects 
during reaction at intermediate energies since interacting isospin asymmetry nuclear 
matter will experience compressed and expanding phases before forming fragments 
around normal density conditions. More complete isospin fractionation processes 
can be obtained from the analysis of the density dependence of neutron and proton 



chemical potentials as above mentioned that the symmetry potential in the chemical 
potential is a repulsive for the neutron and a attractive for the proton , leading to the 
different chemical potential degradations with nuclear density p for the neutron and 
proton in the low density region below normal nuclear density and in the expanding 
process of the colliding system. However in the nonequilibrium expanding process 
the mass flow is determined by the difference in the local values of the chemical po- 
tential and directed from the higher chemical potential region to lower values until 
equilibzation. Above transport effects induce the unstable isospin asymmetric nu- 
clear matter to separate into a neutron-rich low density phase and a neutron-poor 
high density one. Especially quite different chemical potential degradations with 
nuclear density for neutrons and protons below normal nuclear density, induce the 
isospin fractionation [30]. 

In summary, from the calculation results by using IQMD we can get following con- 
clusions: (1) Nn depends sensitively on the isospin effect of two-body coUision and 
weakly on the symmetry potential for the neutron-poor system in the higher en- 
ergy region . (2) In particular, the medium correction of the miclcon nucleon cross 
section enhances the dependence of Nn on the isospin effect of two-body collision 
. (3) MDI produces also an important role for enhancing the dependence of Nn on 
the isospin effect of two-body collision in the medium. (4) On the contrary ,the 
degree of isospin fractionation < {N/Z)gas >b / < {N/Z)uq >b depends sensitively 
on the symmetry potential and weakly on the isospin effect of two-body collision for 
neutron-rich system at lower beam energy region. 
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Figure captions 



Fig.l The nucleoli emission number Nn as a function of the beam energy for JJl^"^ , 
f/s^"" and t/o^™ with a''" and a""*^" for systems ^^Kr+^^Kr and i^25„+4°C„(see 
text). 

Fig.2 The time evolution of the ratio of nuclear density to normal density ^^"^JJ'"^ 
for four different cases for system "^^Kr Kr at E=150 MeV/nucleon and 
b= 4.0 fm (see text). 

Fig.3 The time evolution of the Nn with MDI and NOMDI for four cases (see text) 
and systems ^^Kr Kr at E=100 MeV/nucleon and b= 4.0 fm. 

Fig.4 Time evolutions of an impact parameter average values of the degree of the 
isospin fractionation < {N/Z)gas >b /(< {N/Z)iig >5 (in top window), < 
{N/Z)gas >b and < {N/Z)iiq >b (in bottom window) for the system ^'^^Sn+^'^'^ 
Sn at the beam energy of 50 MeV/nucleon and an impact parameter of 4.0 
fm for the different symmetry potentials C/*^'" and U^^"^ with (t'*° and cr""**" . 



